Coherently enhanced measurements in classical mechanics 
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We show that the recently discovered quantum-enhanced measurement protocol of coherent av- 
eraging [l| that is capable of achieving Heisenberg-limited sensitivity without using entanglement, 
has a classical analogue. The classical protocol uses TV harmonic oscillators coupled to a central os- 
cillator and one measures the signal from the latter. We propose an application to the measurement 
of very weak interactions, and, in particular, a novel route to measuring the gravitational constant 
with enhanced precision. 



A common practice for increasing the signal to noise 
ratio in a measurement is to measure identically prepared 
systems N times and average the measurement results. 
This typically leads to a scaling of the sensitivity (i.e. the 
smallest resolvable change in a parameter that we want 
to measure) as 1/y/N, a scaling that is known in quan- 
tum measurement theory as the "shot noise limit" or the 
"standard quantum limit" (SQL), even though there is 
nothing genuinely quantum about this scaling: It holds 
whenever the central limit theorem applies. 

With the rise of quantum information theory, the 
exciting possibility that the 1/y/N behavior might be 
improved upon has received large attention [2-23]. 
It was shown that if one puts the N systems into an 
entangled state, a scaling as 1/N can be achieved, known 
as the "Heisenberg limit" (HL) Q. Examples include the 
use of NOON states in a Mach-Zehndcr interferometer 



24 l25j , or squeezed spin states for magnetometers based 
on atomic vapors 26]. Unfortunately, the entangled 
states required in these schemes are very unstable and 
prone to decoherence. Experiments with NOON states 
showing a slight improvement over the SQL have not 
surpassed yet the stage of more than a few entangled 
photons [3, [I?]]. In fact, the situation here is much more 
unfavorable than even for a quantum computer: whereas 
for the latter it should be enough to fully control a few 
hundred to a few thousand qubits in order to outperform 
existing classical computers for specific tasks such as 
factoring [28| or data fitting (29|, classical experiments 
such as LIGO have already sensitivities of the order 
10~ 22 /\/Hz [3JJ, |31|. To compete with such performance 
using an etangled state of N particles, one would have 
to entangle a macroscopic number of particles (or create 
a NOON state with a macroscopic number of photons), 
which seems out of reach considering the experimental 
difficulties of creating a NOON state with just 4 photons. 
Also, from theoretical grounds, it has become clear that 
for NOON states the slightest amount of decoherence 
leads back to the SQL scaling for sufficiently large 
N 



The commonly held believe that entanglement is nec- 
essary for reaching the HL is based on propagation of the 
quantum mechanical state of N distinguishable particles 
with a very simple hamiltonian: H = Y^j—i hi(x), where 
hi(x) is a single particle hamiltonian [3j. No interactions 
between different particles are considered. Recently it 
was shown that a system with fc-body interactions offers 
a scaling of the sensitivity as A| J cx 1/N k ^ 1 ^ 2 with- 
out initial entanglement (and 1 /N k with initial entangle- 
ment 0-[li|), even though interactions themselves will 
ultimately have to scale down with TV if the total energy 
is to remain an extensive quantity. The A:-body interac- 
tions typically lead to squeezed states and ressemble in 
this respect the earliest examples of quantum-enhanced 
measurements that proposed the use of squeezed light 
[35L With indistinguishable particles, as obtained 

naturally e.g. from a Bose-Einstein condensate, one may 
avoid entangling the particles as well for surpassing the 
SQL limit, even though the definition of entanglement is 
more tricky in this case 37 1. 

Another possibility is to have N distinguishable 
systems interact with a N + 1st system and read out the 
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-34j. For niche applications, such in biological 
systems that require low intensities, these methods may 
nevertheless be interesting [2fj| . 



This method has the advantage that the 
total system needs to accomodate only N interaction 
terms, such that at least in principle the interaction itself 
may be independent of N. Furthermore, the scaling 
with N appears to be stable under local decoherence, 
and even decoherence itself can be used as a signal, 
if the N + 1st system is an environment. The effect 
can be understood as "coherent averaging" : a phase 
accumulates in the state of the N + 1st system from the 
interaction with the N other systems. No entanglement 
is needed. These properties make one wonder whether 
a classical analogue of this mechanism exists. In the 
present Letter we show that this indeed true. More 
specifically, we show that there is a phase accumulation 
mechanism in the classical motion of a central harmonic 
oscillator interacting with N other harmonic oscillators 
that is completely analogous to the quantum mechanical 
scenario. The phase accumulation allows one to achieve 
a sensitivity that scales as 1/N, just as in the quantum 
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case with Heisenberg-limited sensitivity, even though 
there is of course nothing quantum. Thus, while the 
found sensitivity is, contrary to the quantum case, not 
the expression of a generalized Heisenberg uncertainty 
relation, it shows nevertheless that even in the classical 
realm there are situations where one can improve upon 
the venerated averaging of N independent measurement 
results by coupling the same resources "coherently" to a 
N + 1st system, and measuring the latter. 

Model. Consider a classical harmonic oscillator with 
frequency u>q = harmonically coupled to N other har- 
monic oscillators with frequency uii, i = 1,...,N. The 
Hamilton function (with masses m,- = 1) reads 
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where pi and qi are the canonical momenta and coordi- 
nates, respectively, and £ 2 denotes the coupling strength, 
such that £ has the dimension of a frequency. This is the 
parameter we want to determine. The oscillators need 
not be mechanical oscillators, of course. 

The total potential energy in the problem can be 
rewritten as a quadratic form, 
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The problem can be solved by diagonalizing C . We are 
interested in very small couplings, £ 2 <C w| , i — 0, . . . , N, 
where diagonalization can be performed perturbatively, 
starting from the uncoupled (squared) eigenfrequencies 
A, = uif, i = 0,...,N. Furthermore, we will restrict 
ourselves to the situation where the are narrowly dis- 
tributed about a central frequency uJ, and sufficiently off- 
resonant from f2, such that non-degenerate perturbation 
theory (PT) will suffice to obtain the correction to O. 
To order 0(£ 2 ) we have \ = n 2 + N£ 2 + 0(£ 4 ) and 
A, = cof + £ 2 + O(^), I = 1, . . . ,N. The perturbed 
eigenmodes u; are summarized in the orthogonal trans- 
formation matrix U = (uo, Uj. . . . , Ujy) to order 0(£ 2 ) 
as 
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We consider two sources of uncertainty: i.) an uncer- 
tainty in the frequencies Wi, i = 1, . . . , N, and ii.) time 
dependent noise. 



The equations of mo- 
CiiOo = are de- 



i.) Uncertainty in frequencies. 
tion without external driving, cji + J2j ^ijlj 
coupled into N + 1 independent harmonic oscillators by 
the transformation q > q = U'q. Solving them and 
transforming back leads to the response of the central 
oscillator. In order to simplify expressions, we specialize 
to vanishing initial speeds for all oscillators, <jj(0) = 0, 
j = 0, . ,.,N. We will furthermore assume N£ 2 <C O 2 , 
such that v 7 ^ = n(l + N£ 2 /(2n 2 ))+0(Z 4 )). This limits 
N, but for very small £ 2 , N can become very large (see 
also the comments below for the validity beyond PT). To 
0(£ 2 ) we then get 
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The appearance of a phase shift that scales proportional 
to N and the parameter to be measured is reminiscent 
of phase superresolution [27j j . This signal can be recov- 
ered by mixing the response of the central oscillator with 
a cos(fit) signal corresponding to the unperturbed oscil- 
lator, i.e. one multiplies qo(t) with cos(f2i), which cre- 
ates two signals, one that oscillates with the sum of the 
two frequencies, the other with the difference. The latter 
varies very slowly, and can be isolated with a low-pass fil- 
ter. If we assume the spectral width of the low-pass filter 
to be much smaller than 2f2 and ujj — Q, the remaining 
signal s(t) reads 

s(t) 

where r({uii}) 
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is a random variable whose distribution is given by the 
distribution P({ujj}) of the Wj. The smallest uncertainty 
with which £ 2 can be measured is given by 
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where (. . .) means average over P({uj}), a(s(t)) is the 
standard deviation of s(t) with respect to this distribu- 
tion, and M is the number of measurements. It has the 
meaning of the smallest variation in £ 2 that moves the 
average of the signal at least a distance given by the 
width of the distribution of the signal. A short calcula- 
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tion yields 
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If one waits long enough (N^ 2 t/(2fl) ^> 1), the first 
term in the denominator dominates. If we set in addition 
qo(0) — 0, we obtain the final result 



1 2n 

N J~Mt 



cot 



Nt?t 
2fl 



(r) 



(9) 



The prefactor 1 /N identifies this minimal uncertainty un- 
der the noise process considered as "Heisenberg-limited" . 
Of course, this has nothing to do with Heisenberg's un- 
certainty relation. Rather, we have considered a classical 
noise process (uncertainty in the original frequencies w,), 
but have found a way to reduce the resulting smallest un- 
certainty with which the parameter £ 2 can be measured 
from a 1/yN scaling (that would be obtained by measur- 
ing it separately from each system i and the central oscil- 
lator, and then averaging) to a 1/N scaling. The process 
how this happens is completely analogous to the quan- 
tum mechanical collective phase accumulation described 
in [3^, Hj| : N systems interact with a common central 
system, and lead to an accumulated phase proportional 
to N and the parameter to be measured. This manifests 
itself in an oscillation with a frequency proportional to 
A£ 2 in a corotating frame — the equivalent of homo- 
dyne detection in the quantum optical setting. It also 
leads to the same scaling with t, namely as and not 
the usual l/y/t. It means that the sensitivity per square 
root of Hertz, S^ in y/i, still decreases as l/y/t, just as in 
the Heisenberg limited quantum case. Different from the 
quantum-mechanical case is, however, the restriction to 
N such that iV£ 2 <C SI 2 . This is true beyond the validity 
of PT, as is seen by analysing the very strong coupling 
limit, where C can again be diagonalized analytically. 
Nevertheless, the method proposed here may be advata- 
geous for very small interactions, where N can be very 
large. 

ii.) Time- dependent noise In the presence of time- 
dependent noise forces fi(t) acting on oscillator i, the 
equations of motion in the original oscillator coordinates 
qi read 



9i+^2 C ijQj =/«(*)• 
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After transformation to the eigenmodes q~i defined 
through qj = J2i Ujiqi we get 
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A special solution of this equation can be found with the 
help of the Greens-function of the harmonic oscillator. 



The back transformation to the original coordinates gives 
for the central oscillator 

qo (0) 



qa(t) = go(0)cos(VA i) H 7=sin(vA i) 
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The noise enters here ahead at order £°, i.e. perturbs 
even the uncoupled central oscillator. Nevertheless, we 
will now see that the phase accumulation of the central 
oscillator due to the coupling to the N other oscillators 
still leads to a 1 /N scaling of the sensitivity. 

We restrict ourselves again to qj(0) = for all j 

Y and A£ 2 < fi 2 such that y^O = ^(1 + 

A£ 2 /(2f2 2 ))+0(£ 4 )), and calculate the direct response to 
the noise. For simplicity we consider noise with zero av- 
erage, (fo(t)) = V£, where (. . .) means now average over 
the noise-process. We then have (qo(t)) = qo (0) cos y/Xot 
and a 2 (qo(t)) = a 2 (n(t)), where 



n(t) 



is the noise response. 
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ii.a White noise Consider first white noise, defined 
through (/o(*i)/o(*2)) = foTS(h - t 2 ), where we have 
introduced a unit of time T for dimensional grounds, in 
addition to the force amplitudes fo- One then immedi- 
ately gets 



1 sin 2 



%(t- 



—f-dt' < — 
Ao 
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In fact, for large times, t ^> l/yAo, one has a 2 (qo(t)) ~ 

All the dependence on N of the sensitivity arises 
again from the derivative of (qo(t)} and thus Ao with re- 
spect to £ 2 . Inserting everything in eq.([7j), we are led 
to 
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where for large times still a factor l/y/2 can be gained 
on the rhs. Again, we see that the result scales as 
1/N. The time dependence is different from the 
previous case (and the typical quantum situation at 
the HL): for large times the smallest resolvable <5£ 2 de- 
cays only as 1/y/i, just as in the standard quantum limit. 

ii.b Colored noise. The above considerations are easily 
generalized to colored noise. In fact, unless the noise fo 
on oscillator depends already at order £° on N (which 
appears to be a highly artificial situation, since with- 
out interaction the central oscillator should not "know" 
about the number of additional oscilltors), a 2 (n{t)) is in- 
dependent of N, and the same scaling analysis concerning 
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TV therefore applies and always leads to a 1 /N scaling of 
the sensitivity. Only the time dependence will differ. As 
a more general example, consider stationary colored noise 
with a correlation function (/o(*i)/o(*2)) = /oC(*i — t%) 
where we take C(0) = 1, and C(— t) = C(t). One then 
easily finds the upper bound 
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If the correlation function vanishes for t > t c , one has 
0- 2 (<?o(*)) < U^-b{t) with 
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where we have used \C(t)\ < |C(0)|. Correspondingly, we 
have for the sensitivity 
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which again scales as 1/N. 

A possible application of this "coherent averaging" 
technique might be a novel way of measuring the 
gravitationl constant G, which is one of the least well- 
determined natural constants with a relative uncertainty 
of order 1(T 4 and the 1986 CODATA recommended 
value based on conflicting experimental results [U] . One 
of the reasons for this dire situation is the extremely 
weak strength of the gravitational interaction. This, 
and the impossibility to shield the gravitational field 
from other disturbing bodies, render the determination 
of the absolute value of G very difficult, in spite of 
continued strong interest, driven in part by attempts to 
detect a variation of G as function of distance, time, or 
other physical quantities. Since Cavendish's pioneering 
work in 1798, essentially all lab-experiments attempting 
to measure G were based on a beam balance or a 
torsion pendulum, and measured either a static response 
(some by counterbalancing deflections of the small test 
masses), or a dynamic one (allowing frequency-specific 
analysis synchronized with a periodic excitation, see 
e.g. the recent attempt to measure deviations of the 
1/r 2 behavior below the dark energy length scale of 
about 85/zm Our "coherent averaging" method 

suggests a new, massively parallel way of attempting 
to measure G more precisely: Couple N ^> 1 torsion 
balances graviationally to one central one (consisting 
of N further beams with masses fixed at the ends, and 
all beams attached rigidly to a common axis). Then 
measure the shift in frequency of that central oscillator 
as function of N and the positions of the test masses. 
This will enable a reduction of the uncertainty in the 
measured value of G with a scaling 1/N and should 
become competitive with traditional measurements for 
large N. 



In summary, the above analysis shows that even in the 
classical realm there are situations where it is possible to 
beat the venerated procedure of averaging N measure- 
ment results of a physical quantity that leads to a 1 / \/N 
scaling of the sensitivity. It can be improved upon by 
coupling the N samples "coherently" to a central oscil- 
lator that will pick up a collective phase proportional 
to N and the parameter to be measured, and yield in 
the end a 1/N scaling with the number of samples avail- 
able. Applications might be found in the measurement 
of very weak interactions, such as the gravitational inter- 
action between lab-scale test masses, suggesting a new, 
massively parallel way of determining the gravitational 
constant. 
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